Accurate calculation of polarization-related quantities in semiconductors 
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We demonstrate that polarization-related quantities in semiconductors can be predicted accu- 
rately from first-principles calculations using the appropriate approach to the problem, the Berry- 
phase polarization theory. For IIFV nitrides, our test case, we find polarizations, polarization dif- 
ferences between nitride pairs, and piezoelectric constants quite close to their previously established 
values. Refined data are nevertheless provided for all the relevant quantities. 



The importance of spontaneous and piezoelectric po- 
larization in nitride semiconductors is by now widely 
recognized, and numerous reports of observations and 
practical exploitations of polarization effects in optoelec- 
tronic and. electronic nitride nanostructure devices have 
appearedJj Many experiments have been rather success- 
fully interpreted using recently calculated ab initio val- 
ues of the jSpontaneous polarization and the piezoelectric 
constants.H 

Doubts have recently been raised about the accuracy 
of computed spontaneous polarization values in nitride 
semiconductors. A recent studyu based on a supercell 
methodu reported values of the spontaneous polarization 
of AIN, GaN, and InN that differ wideJy (by -50-100%) 
from those of our recent calculations J3 Our results were 
obtained within density-functional theory using numeri- 
cal methods based on the modern theory of polarization 
and the Berry-phase concept so they are expected to 
be highly reliable. Indeed, they have been used success- 
fully in the experimentalEJ and modelingQ literature. It is 
therefore appropriate to reexamine the accuracy of the 
polarization calculations, to dissipate doubts and confu- 
sion about which values to use, and to clarify the size 
of the expected errors. Since the nitrides provide a se- 
vere test for polarization calculations, we expect our error 
estimates to provide a valid upper bound for all semicon- 
ductors. 

The message of the present paper is fourfold. First, 
we find that it is possible to calculate polarizationa,with 
the accuracy needed to compare with experiment ,£10 pro- 
vided a well-controlled and accurate method is used to 
calculate the polarization, and that the structure is ac- 
curately optimized by state-of-the-art ab initio density- 
functional calculations. The polarization is found to 
be most sensitive to internal structural parameters (in 
wurtzite, the parameter u). Second, absolute polariza- 
tion values resulting from the present revision deviate 
only slightly (^^10%) from the previous values. Moreover, 
observable quantities, namely polarization differences be- 
tween various compounds, are essentially unchanged with 
respect to our previous report .0 Third, we report revised 
values of the piezoelectric constants, again rather close 
to our previous data; we consider that the present values 
are more refined and should be preferred in modeling. 



We accompany the revised piezoconstants by accurately 
calculated elastic constants, providing a self-consistent 
set of data to estimate piezoelectric fields in strained ni- 
trides. Fourth, we briefly address the possible origin of 
the discrepancy between the report of Ref. ^ and our 
results. We suggest that it may be connected with tech- 
nical shortcomings, especially imperfect convergence of 
the equilibrium structure, in the superlattice calculation. 

We analyze the ]D£)larization as obtained using the 
Berry-phase methodH^Ll within two different density- 
functional theory (DFT) exchange-correlation schemes. 
Specifically, as do the authors of Ref. ^, we use the VASP 
packageQ and the pseudopotentials provided therewith. 
We carry out calculations using both the generalized gra- 
dient approximation (GGA) to density-functional theory 
in the Perdew-Wang PW91 version, and the local-density 
approximation (LDA) in the Ceperley-Alder-Pepdew- 
Zunger form (used in Ref. |[). Ultrasoft potentialaa are 
used (Ga and In d electrons are treated as valence) at a 
conservative cut-off of 325 eV, and reciprocal-space sum- 
mation is done on an (888) Monkhorst-Pack mesh. 

Our results are as follows. In Tables |, |l| and III we 
list the structural parameters a, c/a, and ei = u — Uidcai, 
and the spontaneous polarization Pgp, for the III-V ni- 
trides. To obtain them, we optimize the structure within 
both the GGA and LDA, then calculate the sponta- 
neous polarization with the Berry-phase technique. In 
addition, we calculate the polarization for sets of struc- 
tural parameters reported by others. Specifically, we 
use the LDA-pseudopotential lattice paranieters reported 
by Bechstedt, GroBner, and Furthmiillercl (BGF), and 
those calculated with the FLAPW method within the 
LDA by Wei and Zunger (WZ).y The results are labeled 
'Present (LDA)', 'Present (GGA)', 'BGF (LDA)', and 
'WZ (LDA)' in Tables | to g. Next, in Table we give 
the spontaneous polarization differences for some rele- 
vant cases. Finally, in Table ^ we report the piezoelectric 
constants, dynamical Born charges, spontaneous polar- 
izations, and elastic constants (relevant to the piezofields 
in strained nitrides layers) as calculated in the GGA ap- 
proximation, which are proposed as revised values to be 
used in modeling of experiments involving macroscopic 
polarization effects. 

From Tables | to [II , it is evident that the Berry-phase 
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results are quite homogeneous despite the methodological 
differences of the various methods and DFT parameteri- 
zations (the main exception, GaN in the Bechstedt struc- 
ture, is discussed below). As pointed out helow, they dif- 
fer considerably from the values calculatedu in the super- 
cell approach. The major source of these deviations are 
inaccuracies in the determination of the internal param- 
eter u. In the ideal wurtzite structure, w=Uidoai=0.375 
in units of c/a. Here we compare different calculations 
through the parameter ei = u — Uideai, expressed in units 
of 10~^ c/a. It can be noticed from the Tables that po- 
larization is very sensitive to ei, though relatively in- 
sensitive tB|r©ther structural parameters, as also noted 
previouslyB'tll For GaN, the ei value of Ref. ^ deviates 
significantly from the others, and indeed we calculate a 
corresponding Berry-phase Pgp which is more than a fac- 
tor of two larger than all the others. The polarization, of 
— 0.074 C/m^ obtained in the superlattice calculations is 
also too large by a factor 2. This suggests that the origin 
of the discrepancy for GaN is the overestimation of the 
ei paraniater {ei=6.5) in Ref. |[ Since the same group 
reportedli3 ei=1.5 for GaN (in good agreement with other 
entries in our Table and with experiment) using the same 
code and pseudopotentials, but a much finer k-point sum- 
mation mesh, the deviation is possibly due to k-space 
summation. 

For AIN, the largest discrepancy (between the val- 
ues for the GGA and LDA-Bechstedt structures) is 
13mC/m^, i.e., a difference of ^1.0 in ei. Indeed, the 
GGA structure is closer to experiment, and it should be 
more reliable. In any event, the superlattice value of 
-0.120 C/m2 is about 20-25 mC/m^, or 25%, larger than 
all the Berry-phase values. For InN, ei is quite homoge- 
neous throughout, and indeed so are the polarizations, all 
within 5%. The superlattice polarization is — 0.050 C/m^, 
also in the vicinity of the Berry-phase values. 

From the data in Tables ||, III , and |l[ it appears that 
the Berry-phase polarization values are quite stable and 
largely independent of exchange-correlation scheme, once 
the structure (mostly, the ei parameter) is deternuned 
correctly. The polarizations previously determinedcl for 
AIN, GaN, and InN were -0.081, -0.029, -0.032; the 
present GGA values are 11%, 9% and 23% larger. We 
also remark that the GGA values obtained in this pa- 
per produce polarization differences between nitride pairs 
that compare quite well with those in current use.0 Table 
IV reports the comparison. 

We also recalculated the dynamical charges and piezo- 
electric constants of the nitrides in the GGA approxi- 
mation to verify our earlier results. Also, we evaluated 
within the GGA the elastic constants needed in the de- 
termination of strain in epitaxial nitride layers, which 
bears on their built-in piezoelectric fields. The results 
are fairly close to our previously reported values (Table 
II of Ref. ^ . The data listed in Table ^ are proposed as 
updated reference values. [— , 

It has recently been shownEJ that the proper piezoelec- 
tric response is related to a current flow across a sample 



in response to polarization. This response corresponds 
to a set of so-called proper piezoelectric constants. These 
differ from the polarization derivatives with respect to 
strain (the improper constants) by a quantity propor- 
tional to the spontaneous polarization, hence of order 
10%. The piezoelectric constants—just reported in Ta- 
ble are the improper constantsJlj These constants are 
those to be used in the interpretation or modeling of ex- 
periments involving depolpdzing fields and polarization- 
induced interface j&harges.lij On the other hand, in view 
of their definition,ll3 the proper constants are those to be 
compared with experiments based on measurements of 
current across piezoelectric samples aad with values cal- 
culated from linear-response theoryt^ The proper con- 
stants as calculated within the GGA are reported in Ta- 
ble Only the 631 's are affected, while the 633 's are 
unchanged. 

Going back to the accuracy of calculated polarizations, 
it is appropriate to note that the theoretical results are 
relatively stable in other systems also. For BeO (note 
that Ref. ^ misquotes Refs. |^ and ^ about this ma- 
terial), a. longitudinal polarization Pi^— — 15mC/ni'^ is 
reporteco for a length-optimized clamped-ion wurtzite 
half cell. This corresponds to a transverse Pt = 
£00 Pl= — 48mC/rtn^, in excellent agreement with our 
Berry-phase valuaa of — 50mC/m^. Also, the trans- 
verse spontaneous polarization of ZuO was found to 
be — 57mC/m^ p^ing pseudopotentialsB and — 50mC/m^ 
using FLAPW.EiJ Again, the accuracy depends most 
critically on the internal parameters. A similar but. 
more stringent test for GaN/ AIN strained superlatticestll 
showed that the calculated interface charge agrees to 
within less than 1 mC /m^ with that predicted from bulk 
transverse polarizations and dielectric constants.Ej 

We briefly note that the structural issues discussed 
above do not account for all the discrepancies observed 
with the calculations of Ref. ^. Two possible additional 
sources of error are related to the presence of large de- 
polarizing electrostatic flelds in the supercell (see Fig. 3 
of Ref. §). First, the length of the clamped-ion superlat- 
tice was not optimized. In a non-zero field E^, this has 
the same effectEj of an inverse-piezoelectric distortion. 



leading to a spurious polarization SP^ 
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(with 633'' the clamped-ion piezoconstant and C33 the 
elastic modulus). However, this is only a small correc- 
tion (~1 mC/m^). Second, the field causes a large po- 
tential drop over the superlattice, amounting respectively 
to about 1.4, 3.3 and 6.0 eV in InN, GaN, and AIN (in 
the wurtzite region; see Fig. 3 of Ref. ||). This drop 
is larger than the calculated LDA gap of the mateaals. 
While there seems to be no evidence for instability,cl we 
presume this could be a significant source of error. 

In summary we have shown that the polarization of ni- 
tride semiconductors can be calculated within the DFT- 
based first-principles Berry-phase approach in a stable 
and reproducible manner. Calculated polarizations can 
therefore be used safely in comparisons with experiment 
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in nitride quantum structures. Given the demanding test 
case provided by the nitrides, we beheve this conclusion 
apphes to any semiconductor whose structure, and in 
particular whose internal parameters, can be accurately 
predicted. 
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TABLE I. Structure and polarization of GaN. The 
Berry-phase spontaneous polarization Psp has been calcu- 
lated for each structure and setting listed. The structures 
are those otptained in the LDA by Bechstedt, Grofiner, and 
FurthmiilleJa (BGF) and by Wei and Zungeri (WZ), and 
by ourselves (Present) in both the LDA and GGA. Lattice 
constant a in A, ei in 10~^c/a, and Psp in C/m^. Experi- 
mental values from Ref. 0. Ref. I reports a supercell-based 
Psp = -0.074C/m2 for the BGF (LDA) case. 



Structure from 


a 


c/a 


ei 


p 

sp 


BGF (LDA) 


3.f50 


1.6310 


6.5 


-0.080 


WZ (LDA) 


3.f80 


1.6259 


1.8 


-0.032 


Present (LDA) 


3.131 


1.6301 


1.6 


-0.032 


Present (GGA) 


3.197 


1.6257 


1.9 


-0.034 


Experiment 


3.19 


1.627 


2.0 
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TABLE IL Structure and polarization of AIN. Details as 
in Table |. Ref. |^ reports a supercell-based Psp=— 0.120 C/m^ 
for the BGF (LDA) case. 



Structure from 


a 


c/a 


ei 


p 

sp 


LDA (BGF) 


3.080 


1.6070 


7.4 


-0.103 


LDA (WZ) 


3.112 


1.6009 


6.9 


-0.094 


LDA (present) 


3.070 


1.5997 


7.1 


-0.099 


GGA (present) 


3.108 


1.6033 


6.4 


-0.090 


Experiment 


3.11 


1.601 


7.1 





TABLE in. Structure and polarization of InN. Details as 
in Table |. Ref. |^ reports a supercell-based Psp =— 0.050 C/m^ 
for the BGF (LDA) case. 



Structure from 


a 


c/a 


ei 


p 

sp 


LDA (BGF) 


3.530 


1.6320 


5.0 


-0.043 


LDA (WZ) 


3.544 


1.6134 


4.0 


-0.042 


LDA (present) 


3.509 


1.6175 


3.6 


-0.041 


GGA (present) 


3.580 


1.6180 


3.7 


-0.042 


Experiment 


3.544 


1.613 







TABLE IV. Modulus of the spontaneous-polarization dif- 
ferences (mC/m^) between binary nitrides. The values of 
Ref. ^ and Ref. ^ are reported together with the present GGA 
result. 



1 API 


Present (GGA) 


Ref. 2 


Ref. 3 


AIN /GaN 


55 


52 


46 


GaN/InN 


7 


3 


24 


InN/ AIN 


48 


49 


70 
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TABLE V. Spontaneous polarization (C/m^), improper 
piezoelectric constants (C/m^), dynamical charges, and elas- 
tic constants (GPa) for the nitrides as obtained in the GGA 
approximation. 





p 

^ sp 


633 


631 


Z* 


C33 


C31 


AIN 


-0.090 


1.50 


-0.53 


2.65 


377 


94 


GaN 


-0.034 


0.62 


-0.34 


2.76 


354 


68 


InN 


-0.042 


0.81 


-0.41 


3.10 


205 


70 



TABLE VI. Proper piezoelectric constants (C/m^) for the 
nitrides in the GGA approximation. 





633 


631 


AIN 


1.50 


-0.62 


GaN 


0.62 


-0.37 


IiiN 




-0.45 
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